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ABSTRACT 



Aims. We study the impactor flux and cratering on Pluto and Charon due to the coUisional evolution of Plutinos. Plutinos are those 
trans-Neptunian objects located at ~ 39.5 AU, in the 3:2 mean motion resonance with Neptune. 

Methods. To do this, we develop a statistical code that includes catastrophic collisions and cratering events, and takes into account the 
stability and instability zones of the 3:2 mean motion resonance with Neptune. Moreover, our numerical algorithm proposes different 
initial populations that account for the uncertainty in the size distribution of Plutinos at small sizes. 

Results. Depending on the initial population, our results indicate the following. The number of D > I km Plutinos streaking Pluto 
over 3.5 Gyr is between 1 271 and 5 552. For Charon, the number of D > I km Plutino impactors is between 354 and 1 545. The 
number of D > I km craters on Pluto produced by Plutinos during the last 3.5 Gyr is between 43 076 and 113 879. For Charon, the 
number of D > I km craters is between 20 35 1 and 50 688. On the other hand, the largest Plutino impactor onto Pluto has a diameter 
between ~ 17 and 23 km, which produces a crater with a diameter of ~ 31 - 39 km. In the same way, the largest Plutino impactor 
onto Charon has a diameter between ~ 10 and 15 km, which produces a crater with a diameter of ~ 24 - 33 km. Finally, we test the 
dependence of results on the number of Pluto-sized objects in the Plutino population. If 2 Pluto-sized objects are assumed in the 3:2 
Neptune resonance, the total number of Plutino impactors both onto Pluto as Charon with diameters D > 1 km is a factor of ~ 1.6 - 
1.8 larger than that obtained considering only 1 Pluto-sized object in this resonant region. 

Conclusions. Given the structure of the trans-Neptunian region, with its dynamically different populations, it is necessary to study 
in detail the contribution of all the potential sources of impactors on the Pluto-Charon system, to obtain the main contributor and the 
whole production of craters. Then, we will be able to contrast those studies with observations which will help us to understand the 
geological processes and history of the surface of those worlds. 

Key words, methods: numerical - Kuiper Belt: general 
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1. Introduction 

Pluto and Charon are members of a vast population of icy bod- 
ies beyond Neptune and constitute in fact the first discovered 
binary trans-Neptunian object (TNO). In the trans-Neptunian re- 
gion there are at first four dynamical classes (Chiang et al. 2007 ). 
The Classical Objects those with semimajor axis a greater than 
~42 AU and low eccentricity orbits, the Scattered Disk Objects 
(SDOs) with perihelion distance q > 30 AU with large eccen- 
tricities, the Resonant Objects those in mean motion resonance 
with Neptune and the Centaurs with perihelion distance ^ < 30 
AU. 

The 3:2 mean motion resonance with Neptune, located at 
~ 39.5 AU, is the most densely populated one in the trans- 
Neptunian region. The residents of this resonant region are usu- 
ally called Plutinos because of the analogy of their orbits with 
that of Pluto, which is its most representative member Aside 
from Pluto and its largest moon Charon, the Minor Planet Center 
(MPC) database contains ~ 200 Plutino candidates. 

Pluto and Charon have been exposed to impacts with minor 
bodies as all the objects in our Solar System. Cratering is one 
of the most important processes that determine the morphology 
of the surface of a Solar System object. The understanding and 
quantification of the impactor source population onto an object 
and the observation of the object surface help to understand the 
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dynamical and physical history of both the impactor population 
and the target. 

All of the detailed knowledge of the surface composition of 
Pluto and Charon has been obtained from telescopic observa- 
tions of the spectrum of sunlight reflected from their surfaces. 
Pluto's reflectance spectrum shows absorption bands of methane 
ice, and an absorption band that could be related to the presence 
of CO and nitriles. The surface of Charon can be modeled by 
pure water ice darkened by a spectrally neutral continuum ab- 
sorber (Protopapa et al. l2008l l. 

However, it will not be until 2015 that we will have a real 
idea of the morphology of the surfaces of Pluto and Charon, with 
the fly by of the Nasa's New Horizon Pluto-Kuiper Belt mission 
to Pluto system. The New Horizons mission is the first one to 
the Pluto system and the Kuiper Belt and was launched on 19 
January 2006 on a Jupiter Gravity Assist trajectory toward the 
Pluto system for a 14 July 2015 closest approach. It will study 
the Pluto system over a 5-month period beginning in early 2015 
in particular providing measurement of cratering records. 

The comparison of the predicted theoretical crater produc- 
tion from a given source with the observed surface of Pluto and 
Charon may account for the geological processes acting on the 
surface of the objects and if cratering collisions onto Pluto and 
Charon are an important surface modification process. 

Then, it is very important to study all the possible sources of 
crater production on Pluto and Charon in order to know the total 
crater production and to contrast them with observations. 
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The possible main contributors to the impactor flux on Pluto 
and Charon would be on one hand all type of comets. This topic 
has been addressed by Weissman & Stern (1994), Durda & Stern 
( I2OOOI 1. and Zahnle et al. ( I2003I I. In general comets have eccen- 
tric orbits that can cross Pluto's orbit when they enter to the plan- 
etary region from their source, either the trans-Neptunian zone 
or the Oort Cloud, or when they leave the inner Solar System 
because of the perturbations of the planets. 

Weissman & Stern ( |1994| l estimated current impact rates of 
comets on Pluto and Charon. They showed that cratering on 
both bodies is dominated by Kuiper Belt and inner Oort cloud 
comets. Then, Durda & Stern (2000) calculated collision rates in 
the Kuiper Belt and Centaur region through a numerical model. 
They estimated that the flux of Kuiper Belt projectiles onto Pluto 
and Charon is ~ 3 - 5 times that of Weissman & Stern (119941 ). 
Later, Zahnle et al. (2003) studied the cratering rates for the 
moons of the jovian planets and Pluto produced mainly by eclip- 
tic comets, obtaining results consistent with the previous esti- 
mates. We will analyze those studies in a discussion section. 

On the other hand, Plutinos may be the other important 
source of impactors on Pluto and Charon. Plutinos shared the 
same dynamical conditions that Pluto and Charon, in fact they 
are all located in the same 3:2 mean motion resonance with 
Neptune so they have a certain collision probability. Recently, 
de Elia et al. (^2008^ performed a collisional evolution of Plutinos 
and determined collisional rates among these objects. That work 
allows us to study the impactor flux on Pluto and Charon due 
to the collisional evolution of Plutinos. Then, in this paper we 
are going to evaluate the contribution of Plutinos to the impactor 
flux and cratering on Pluto and Charon and also determine if the 
Plutino population can be considered a primary source of im- 
pactors on the Pluto-Charon system. 
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Fig. 1. Impact strength. The dashed lines represent the two dif- 
ferent Qs laws used in our simulations. The Qd law from Benz 
& Asphaug ( 119991 ) for icy bodies at 3 km s ' is plotted as a solid 
line. 



made by de Elia et al. (120081 ) we use in our model three different 
initial populations, which are defined as follows 

- Initial Population 1, with a cumulative power-law index p of 
3.0 for D< 60 km, 

- Initial Population 2, with a cumulative power-law index p of 
2.7 for D < 60 km, 

- Initial Population 3, with a cumulative power-law index p of 
2.4 for D < 60 km. 



2. The Full Model 

To simulate the collisional and dynamical evolution of the 
Plutino population, we use the statistical code developed by de 
Elia et al. (2008). This algorithm considers catastrophic colli- 
sions and cratering events, and takes into account the main dy- 
namical characteristic associated to the 3:2 mean motion reso- 
nance with Neptune. In the following, we give a brief descrip- 
tion of the initial populations, the collisional parameters and the 
main dynamical considerations used in our model. 

2.1. initial Popuiations 



From Kenyon et al. ( 120081 ). the cumulative size distribution of 
the resonant population of the trans-Neptunian region shows a 
break at a diameter D near 40 - 80 km. Moreover, for larger 
resonant objects, the population seems to have a shallow size 
distribution with a cumulative power-law index of ~ 3. From 
this, the general form of the cumulative initial population used 
in our model to study the collisional and dynamical evolution of 
the Plutinos can be written as follows 



A^(> D) = Ci( ^^^ I forD < 60km, 



/lkm\^^ 

/lkm\^ 
A^(> D) = C2 forD > 60km, 



(1) 



V^-3 



where C2 adopts a value of 7.9 x 10^ and Ci = C2 (60)'' 
by continuity for D = 60 km. The size distribution of Plutinos 
at small sizes is uncertain and one can find different proposed 
power law indexes in the literature. Then, following the analysis 



2.2. Coiiisionai Parameters 

Here, we adopt constant values of the intrinsic collision prob- 
ability (Pic) and the mean impact velocity (V) for Plutinos 
derived by Dell'Oro et al. (2001"). Based on a sample of 46 
Plutinos, these authors computed values of (Pic) and (V) of 4.44 
+ 0.04 X 10"^^ km"- yr"' and 1.44 + 0.71 km s"', respectively. 

As for the impact strength, O'Brien & Greenberg (120051) 
showed that the general shape of the final evolved asteroid pop- 
ulation is determined primarily by the impact energy required 
for dispersal Qo, but variations in the shattering impact specific 
energy Qs and the inelasticity parameter /^^ can affect such a 
final population even if Qd is held the same. According to these 
arguments, we decide to choose a combination of the parameters 
Qs and fke that yield the Qo law from Benz & Asphaug (119991 ) 
for icy bodies at 3 km s"^ 

de Elia & Brunini (120071 ) analyzed recently the dependence 
of their numerical simulations on the shattering impact specific 
energy Qs . According to this work, the smallest gaps between 
Qs and Qo curves lead to the smallest wave amplitudes in the 
size distribution of the final evolved population as well as to the 
highest ejection rates of collisional fragments. Moreover, that 
study also indicates that the formation of families is more efifec- 
tive for the simulations with a small gap between Qs and Qd 
laws. Following these arguments, in this work we decide to use 
two Qs laws, Qs,i and Qs.i, with a small and a large gap with 
respect to the Qd law from Benz & Asphaug ( 119991 ) for icy bod- 
ies at 3 km s ', respectively. The results discussed in this work 
are those obtained using the 2^ 1 law. In Sect. 3.4, we develop 
numerical simulations using the Qso law in order to test the de- 
pendence of our results on such collisional parameter Figure [1] 
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Table 1. Collisional evolution results about Pluto. N(D > 1) 
is the total number of D > 1 km Plutinos impacting Pluto over 
3.5 Gyr, t(D > 1) the time scale on which such impacts occur, 
DuMax the largest Plutino expected to have struck Pluto during 
the past 3.5 Gyr, C(D > 1) the total number of D > 1 km craters 
on Pluto produced by Plutino impacts over 3.5 Gyr, and Z)c,Max 
the largest crater diameter on Pluto. 





Initial Pop. 1 


Initial Pop. 2 


Initial Pop. 3 


N(D > 1) 


5 552 


2 843 


1271 


t(D > 1) (Myr) 


0.63 


1.23 


2.75 


Oi.Max (km) 


22.5 


20.1 


16.8 


C{D > 1) 


113 879 


76 726 


43 076 


£>c,Max (km) 


38.8 


35.5 


30.1 



shows the two Qs laws used in our simulations as well as the Qo 
law from Benz & Asphaug (,1999j for icy bodies at 3 km s"' . 

Once the Qs law is specified, we fit the ineslaticity parameter 
fke to get the Benz & Asphaug ( 119991 ) Qd law. According to 
O'Brien & Greenberg (2005) , we express the parameter /^^j, as 



./k 



ke 



■^'''° 1 1000 km) 



where /keo is the value of /ke at 1 000 km and y is a given ex- 
ponent. Our simulations indicate that the Qd law from Benz & 
Asphaug (I1999i) for icy bodies at 3 km s"' is obtained with good 
accuracy from the combination of the selected Qs law and /^e, 
with /ke„ - 0.27 and -y - 0.7. Such values are consistent with 
those from Davis et al. (I1989I) . 

2.3. Dynamical Considerations 

To study the orbital space occupied by the Plutino population, 
we develop a numerical integration of 197 Plutino candidates 
extracted from the Minor Planet Center database with semima- 
jor axes between 39 and 40 AU. These objects are assumed to 
be massless particles subject to the gravitational field of the Sun 
(including the masses of the terrestrial planets) and the perturba- 
tions of the four giant planets. The simulation is performed with 
the simplectic code EVORB from Fernandez et al. (2002)). The 
evolution of the test particles is followed for 10^ years which is 
a timescale greater than any secular period found in this reso- 
nance (Morbidelli 119971 ). From this, we build maps of the dis- 
tribution of Plutinos in the orbital element planes (a,e) and (a,/), 
which allows us to determine the main stability regions of the 
3:2 Neptune resonance. Such maps are used to assign a charac- 
teristic orbit for every colliding Plutino and to specify the final 
fates of the different fragments generated in the collisional evo- 
lution. A detailed description of this procedure can be found in 
de Elia et al. (i2008) . 



3. Results 

The previously described collisional code allows us to calculate 
the collisional rates of Plutinos onto Pluto and Charon. From 
this, it is possible to calculate the impactor flux of Plutinos of 
different sizes on the Pluto-Charon system. Besides, using a suit- 
able expression we can calculate the crater diameters produced 
by Plutinos on Pluto and Charon. We will present here our main 
results concerning the impactor flux and cratering onto Pluto and 
Charon due to the collisional evolution of the Plutino population. 
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Fig. 2. Cumulative number of Plutino impactors onto Pluto a) 
and Charon b) over the last 3.5 Gyr as a function of impactor 
diameter 

Table 2. Collisional evolution results about Charon. N(D > 1) 
is the total number of D > 1 km Plutinos impacting Charon over 
3.5 Gyr, t{D > 1) the time scale on which such impacts occur, 
Di^Max the largest Plutino expected to have struck Charon during 
the past 3.5 Gyr, C(D > 1) the total number of D > 1 km craters 
on Charon produced by Plutino impacts over 3.5 Gyr, and D^Max 
the largest crater diameter on Charon. 





Initial Pop. 1 


Initial Pop. 2 


Initial Pop. 3 


N{D > 1) 


1545 


791 


354 


t(0 > 1) (Myr) 


2.27 


4.42 


9.9 


i'l.Max (km) 


14.6 


12.5 


9.8 


C(D > 1) 


50 688 


34 541 


20 351 


i'cMiix (km) 


32.7 


28.9 


24 



3.1. Impactor Flux onto Pluto and Charon 

Figure |2] a) shows the cumulative number of Plutino impacts 
onto Pluto over the last 3.5 Gyr as a function of impactor di- 
ameter, obtained from the three different initial populations de- 
fined in Sect. 2. 1 . Moreover, Table [1] summarizes some of our 
results concerning the impactor flux onto Pluto due to the colli- 
sional evolution of Plutinos. From this, the number of D > I km 
Plutinos striking Pluto over 3.5 Gyr is between 1 271 and 5 552, 
while the largest Plutino expected to have impacted Pluto during 
the past 3.5 Gyr had a diameter of ~ 17 - 23 km, depending on 
the initial size distribution. 
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Table 3. Values of the density p, the surface gravity g, and the 
encounter velocity v, used for Pluto and Charon. 
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umulative number of craters on Pluto a) and Charon b) produced by 
Plutinos over the last 3.5 Gyr as a function of crater diameter. 



On the other hand, Figure|2]b) shows the cumulative number 
of Plutino impacts onto Charon over the last 3.5 Gyr as a func- 
tion of impactor diameter Moreover, results about the impactor 
flux onto Charon due to the collisional evolution of Plutinos are 
summarized in Table |2] From this, the number of D > I km 
Plutinos striking Charon over 3.5 Gyr is between 354 and 1 545, 
while the largest Plutino expected to have impacted Charon dur- 
ing the past 3.5 Gyr had a diameter of ~ 10 - 15 km, depending 
on the initial size distribution. 

3.2. Cratering on Pluto and Charon 

Using the impactor fluxes previously estimated, it is possible to 
obtain the number of craters on Pluto and Charon over the last 
3.5 Gyr as a function of crater diameter. From Holsapple (I1993I I. 
the diameter Dc of a hemispherical crater is given by 



De = 1.26A(Ap,/A)"^(1.61gDi/vf) 



i/V 



„2^-«/3 



(3) 



where Di is the impactor diameter, pi the impactor density, pt 
the target density, g the surface gravity of the target, Vi the en- 
counter velocity, and A and a are constants which depend on the 
mechanical properties of the target material. We adopt the values 
of A = 0.2 and a - 0.6 given by Holsapple ([1993) for water ice. 
The encounter velocity Vi is given by 




where U is the hyperbolic encounter velocity, which is equal to 
the mean impact velocity given in Sect. 2.2, and Vesc is the escape 
velocity of the target. Values of the density, the surface gravity, 
and the encounter velocity used for Pluto and Charon are shown 
in Table |3] 

Figure[3]a) shows the cumulative number of craters on Pluto 
produced by Plutinos over 3.5 Gyr as a function of crater di- 
ameter From this, the number of D > 1 km craters on Pluto is 
between 43 076 and 113 879, while the largest crater diameter is 
between ~ 3 1 and 39 km, depending on the initial size distribu- 
tion. 

At the same way. Figures [3] b) shows the cumulative num- 
ber of craters on Charon produced by Plutinos over 3.5 Gyr as 
a function of crater diameter. From this, the number of D > 1 
km craters on Charon is between 20 351 and 50688, while the 
largest crater has a diameter of ~ 24 - 33 km, depending on the 
initial size distribution. 



3.3. Dependence of results on the number of Pluto-sized 
objects 



Vi = (t/2 + vl 



1/2 



(4) 



Kenyon et al. (2008*) suggested that the resonant trans-Neptunian 
population has ~ 0.01 - 0.05 Me in objects with D > 20 - 40 
km. From these estimates, it is possible to infer the existence of 
5 Pluto-sized objects in the whole resonant population. If these 
5 Pluto-sized objects were all in the 3:2 Neptune resonance, we 
would have an upper limit for the large objects in this resonance. 
Brown (I2008I I. based on the completeness of the current surveys 
argued that two or three more large KBOs are likely awaiting 
discovery. Since the existence of 5 Pluto-sized objects in the 3:2 
Neptune resonance would seem to be an overestimation, we as- 
sume an upper limit of 2 Pluto-sized objects in this resonant re- 
gion in order to analyze the dependence of our simulations on 
the number of large objects in this resonance. 

If 2 Pluto-sized objects are assumed in the 3:2 Neptune res- 
onance, the total number of Plutino impactors both onto Pluto as 
Charon with diameters D > 1 km is a factor of ~ 1.6-1.8 larger 
than that obtained considering only 1 Pluto-sized object in this 
resonant region. 



3.4. Dependence of results on the shattering impact specific 
energy Qs 

To test the dependence of our results on the shattering impact 
specific energy Qs , we carry out numerical simulations using the 
Qs.i and Qso laws, which show a small and a large gap with re- 
spect to the Qd law from Benz & Asphaug (1999 J for icy bodies 
at 3 km s"', respectively. We obtain that the Qj,2 law produces a 
size distribution with a wave amplitude larger than that the one 
obtained using the Qs^\ law as was noticed by de Elia & Brunini 
(120071) . However, the whole impactor flux is almost preserved 
using both laws. 
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4. Discussion 



Durda & Stern ( I2000I I developed a model of collision rates in 
the present-day Kuiper Belt and Centaur region. That is a static, 
multizone and multi-size-bin collision rate model that calculates 
instantaneous collision rates on Kuiper Belt objects, Centaurs 
and the Pluto-Charon system. 

The model of Durda & Stern (I2000I I defines a disk in the 
30-50 AU zone whose surface mass density follows a profile of 
the form E oc R^^, with /? = -2. Moreover, the disk adopts an 
average eccentricity (e) and assumes an equilibrium condition 
where the average inclination (/) - |(e). The values adopted 
for (e) are of 0.0256 and 0.2048, which are representative of 
the values observed in the Kuiper Belt at that time. Once the 
disk properties are specified, the disk is binned into a series of 
radially concentric tori 1 AU in width in each of which the size 
distribution of the population is represented by a two-component 
power law of the form N{D) cc D'^dD, where b - -3 for D < 
10 km and b - -4.5 for larger bodies. Durda & Stern ( 120001 1 
carried out their runs assuming the existence of 7 x 10'*-1.4 x 
10^ objects with a radius r > 50 km between 30 and 50 AU, 
according to the observational evidence from Jewitt et al. di 9981 1 
and Gladman et al. (I1998I I. Moreover, they suggested that about 
40 % of the total population are in or near the 3:2 mean motion 
resonance with Neptune. From this, the flux of r = 1 km Kuiper 
Belt projectiles onto Pluto and Charon over 3.5 Gyr is found to 
be approximately 8.9 x lO-' and 1.1 x 10^, respectively. 

However, the observational advances that have occurred dur- 
ing the last ten years have revealed a more complex dynamical 
structure of the trans-Neptunian region. In particular in the re- 
gion between 30 and 50 AU there are very different dynamical 
classes of TNOs, like the classical, resonant, and scattered-disk 
populations. They show different dynamical features and ranges 
of orbital elements as well defined semimajor axis zones, eccen- 
tricities as high as ~ 0.45 and inclinations as high as ~ 45°. 
Moreover, current data provide clear evidence for diff'erences in 
the mass and size distribution parameters among such dynami- 
cal classes (Kenyon et al. 2008). The results by Durda & Stern 
(I2000I I are obtained for objects in that region, but they don't ac- 
count for the different dynamical populations. Particularly, while 
the Plutino population is included quantitatively in Durda & 
Stern's (2000) model, they did not account for the dynamical 
properties of such resonant population. Then, their results should 
be taken with caution. 



5. Conclusions 

We have presented a study aimed at analyzing the collisional and 
dynamical evolution of the Plutinos. Assuming the existence of 
one Pluto-sized object in the 3:2 Neptune resonance, our main 
results are the following. 

Our results depend strongly on the initial size distribution of 
Plutinos. Depending on the initial population, we have obtained 
that the number of D > 1 km Plutinos streaking Pluto over 3.5 
Gyr is between 1 271 and 5 552. For Charon, the number of D > 
1 km Plutino impactors is between 354 and 1 545. Moreover, the 
number of D > 1 km craters on Pluto produced by Plutinos dur- 
ing the last 3.5 Gyr is between 43 076 and 113 879. For Charon, 
the number of D > 1 km craters is between 20351 and 50688. 
On the other hand, the largest Plutino impactor onto Pluto has a 
diameter between ~ 17 and 23 km, that produces a crater with 
a diameter of ~ 3 1 - 39 km. The largest Plutino impactor onto 
Charon has a diameter between ~ 10 and 15 km, that produces a 
crater with a diameter of ~ 24 - 33 km. 



We have tested the dependence of results on the number of 
Pluto-sized objects in the Plutino population. The number of 
Plutino impactors on Pluto and Charon obtained considering 2 
Pluto-sized objects is a factor of ~ 1.6 - 1.8 larger than that 
obtained considering only 1 Pluto-sized object in this resonant 
region. 

Besides, using two different Qs laws, with a small and a large 
gap with respect to the Qd law, we have obtained almost the 
same impactor flux onto Pluto and Charon. 

The complex structure of the trans-Neptunian region, with 
its dynamically different populations, requires a detailed study 
of the contribution of all the potential sources of impactor on 
the Pluto-Charon system, to obtain the main contributor and the 
whole production of craters. When Nasa's New Horizon Pluto- 
Kuiper Belt mission flies by Pluto system we will have images 
of the surface and cratering records could be measured. Then, 
we will be able to contrast the theoretical studies of the produc- 
tion of craters with the observations of the New Horizon mission 
which will help us to understand the geological processes and 
history of the surfaces of those worlds. In a future work we will 
calculate the contribution of the different dynamical classes of 
the trans-Neptunian region to the flux of projectiles onto Pluto 
and Charon. From this, it would be possible to specify the pri- 
mary source of impactors on such bodies and to determine if the 
Plutino population can be considered a significant source. 
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